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Abstract 

The spectrum of a single hole is calculated within the spin-hole model 
using a variational method. This calculation is done for any rotational in- 
variant magnetic background. We have found that when the magnetic back- 
ground changes from a disordered to a locally ordered state, the spectrum 
changes qualitatively. We have also found that the spin pattern around the 
hole is polarized. This problemis related to the study of copper oxide planes 
C11O2 doped with a small number of holes. 
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1 Introduction 



The study of the motion of one hole in strongly correlated systems is very 
important because, in principle, it may contain some features of the proper- 
ties of the C11O2 planes doped with a small number of holes. The one hole 
motion in a quantum antiferromagnetic (AFM) background has been stud- 
ied by several authors JI], 0, [| [3J] in the framework of the t-J model. These 
studies were carried out within an AFM ordered state. The main result of 
these analysis is the existence of a well-defined quasiparticle description for 
a coherent hole propagation with the band minima at (±7r/2, ±vt/2) and the 
bandwidth W of the order of J and not t. These results are in good agree- 
ment with numerical exact diagonalization for small systems || . All these 
investgations concern the magnetic ordered case, and they were carried out 
in the framework of the t-J model which has the difficulty of the constraint 
of no double occupancy at the same site. In this paper we return to the 
Hamiltonian Hp} d derived in, e.g., Refs. J7|, [5|, [|, |K]]. This Hamiltonian is 
obtained from the Emery model by means of a perturbative expansion up to 
the second order in t pd (the hopping parameter between Cu and O sites). We 

(2) 

point out that the t term of the t-J Hamiltonian is obtained from H^_ d by 
projection on local singlet states^]. We will not consider here the fourth or- 
der term which gives the super-exchange energy between spins located on the 
Cu sites. This term will not be relevant in our calculation since we want to 

study essentially the motion of a hole in a disordered magnetic background. 

(2) 

In the framework of the H p _ d Hamiltonian we will compute the dispersion 
relation E(k) by applying a trial wave function. The results agrees with pre- 
vious calculation for the t-J model when we consider the limit of an ordered 
magnetic background. 

Our starting point is the H^} d Hamiltonian (for shortness we set H as 

(2) 

Hp_ d ) which can be written as 

H =J2 T ij fiafja + 9 J2 W ij W if fja W fj'a'A (1) 
i,j,cr i,j,j' t a,<j' 

where (f^, fi a ) is a fermion field with the spin index a . Si is the spin oper- 
ator on the Cu site i, f = (r x , r y , r z ) are the Pauli matrices and 

T* = T Tr £ £ « eik(i ~ j) ( 2 ) 
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^ = I E "k ^~ 3) (3) 

where Ek = 4(1 — 7^.), wt = yfe^ and 7^ = (cosfe^ + cosk y )/2. The parame- 
ters T and g are related to the parameters A = e p — e^j Ud and t p( j by the 
following relations 

1 . . 1 1 



The Hamiltonian ([!]) describes a system of itinerant holes interacting with 
a magnetic background. We notice that if we restrict the g term to a local 
spin-spin interaction (i.e. take Wij = wo5ij) we get the phenomenological 



Hamiltonian considered by Monthoux and Pines[13 



2 Construction of the trial state for the one 
hole motion 

The most general state for a single hole is 

ij = £ fl A{ia) $ (6) 

where $ is some magnetic background and A(ia) is some combination of 
spin operators. The meaning of Eq. ([]) is clear : the motion of the hole will 
perturbe the magnetic background generating all sorts of excited magnetic 
states provided by the operator A(ia). 

In what follows, we will suppose that the magnetic background <3> is rotational 
invariant,^, e. 

Stot® = 

For the Hamiltonian under investigation Eq.([TJ), if we set g = 0, the following 
state for one hole with spin up 

E i>o(i)fh® (7) 
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is an eigenstate for the T term of the Hamiltonian Eq. QT]) , if if) (i) = e lk ' 1 for 
a given k. On the other hand the g term can be rewritten in the following 
form 
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i,3,3 



(8) 



In this form we have splitted the Hamltonian into two parts, one which 
contains singlet operators and the other one which contains triplet operators. 
These singlet and triplet operators are more general than those considered 
by Zhang and Rice[|l!|, because here they are non local. 
From Eq.(g) we see also that the singlet states and the triplet states are 
separated by an energy 2g. If we want to look for low energy excitations and 
we neglect mixing between the two bands, we can write the following state 
for the one hole with spin up 



(9) 



Then a natural trial state for the full Hamiltonian Eq. (|1|) can be written as 



da® 



(10) 



To simplify calculation we take ipi(i,j) = (pi(J)Wji. Then the state fllCf ) can 
be rewritten in a more transparent way 



11^ 



i [ 3 

where (p D (i) = if) (i) - \ J2j Vi(J) w ij- 

Examples for the magnetic state $ are the singlet Resonating- Valence- 



Bond states considered by Anderson, Dougot and Liang |14H which may be 
long ranged or short ranged ordered. 
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3 Dispersion relation for a single hole 



The dispertion relation E(k) will be calculated by using the variational 
method. We will minimize (ip,Hip) keeping the norm fixed. Straight- 

forward calculation, using the property : J2i w u w ij — T^, leads to 



{ip : H^) = T i:j < (p (i)(po(J) + w u w ji' ViWvi (0 < > \ + 

i,j { 1,1' J 

+g Tu'Wji {(pi(l')<p (j) < SySi> +h.c} - 

-gJ2 T u > fflf)<Pi{l') < SiS v > +h.c} + 
i,v 

i,v 

+ ig T hl3 f hh ^Opi(Z 2 ) < Sl.(Sl x Si) > (12) 

with = Tij/T and the average < ... > defined as 

< S a ...S b >= ($,S a ...S 6 $) (13) 

Now if we suppose that the magnetic background $ [in Eq.flTTDI is invari- 
ant under the time reversal operation K , we will have 

(K$,S* h S v h S\KH>) = (®,S x h S y l2 S z l3 $) 

= -(<Z>,S x h S y l2 S z l3 <Z>) (14) 

for any Z 1; l 2 , I3 (all different), since KSK^ 1 = —S. And finally we have 

< Sl-iSl x 4) >= (15) 

Consequently the last term of Eq. flT2| ) vanishes. 

The energy will be calculated by minimizing the energy functional 

W(V,A) = (^ J ffV)-A(V,V') (16) 
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Using Eq.(|T2|) and 



tyrf) = £ | <p (i) | 2 + E T «'ViO>iO") < ^A' > (17) 



3 J 



we get the following equation for the energy E 
I Te k - E gw k a k 

\ gw k a k Tc k + gd k - Ea } 




(18) 



where (p (k) and <fi(k) are the Fourier components of ip (i) and fi(i) respec- 
tively. On the other hand 

' ak = 3 - 4xoi 7fc ; c fc = 15 - 32xoi Ik + 4x o2 (47| - 1) + ij 2 k(xu ~ X02) 

< 

4 = 12 + Xoi - 8a fc 

(19) 

with 

Xoi = - < S Q S ei > ; xn =< S S ei+e2 > and X02 = - < S S 2ei > (20) 

ei and e 2 denoting the two unit vectors of the square lattice. From Eq.fll8D, 
the lowest energy band is given by 

E k /g = I Utk + -) + -- 
z 1 a k a k 



- J [v(s k --)--? + 4e k a k \ (21) 
V a k a k 

with i] = T/g and for a given k we have the following solutions for cp (i) and 
(p {i;k) = w k a k U(k)e ik - 1 

Vl (i;k) = -( V e k -E k /g)U(k)e ikA (22) 
where the function U(k) is determined by writing (if>, if)) — 1 and we get 

U{k? = e k al + l[ V l k -E k /gf (23) 
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The dispersion relation E k depends only on the the spin-spin correlation func- 
tion for nearest and next nearest neighbors. This is due to our expression for 
the trial state. Introducing more and more spin excitations would essentially 
produce new spin correlation functions at larger distance. Hower since we 
are interested essentially in a background with short correlation functions, 
this would merely not chane our result at least qualitatively. Knowing Xoi, 
Xn and X02 we evaluate then the energy E k . 

If we set rj = 0, we can have a qualitative idea about the positon of the 
minimum of Ek 

E k /9 = -\{\ — I +J\-\ 2 +4e k a k ) (24) 
2 [ a k V a>k J 

It is easy to see that \ d k /a k \ is maximal for k = (0, 0) and Ae k a k is maximal 
for k = (tt,7t). Then since to lower E k both of the quantities should be 
maximal, there will be a competition between them. For Xoi = , \ d k /a k \ = 
4 and consequently the minimum will be at k = (it, it). For Xoi 7^ the 
position of the minimum will be in between (0, 0) and (tt, tt). 

3.1 Paramagnetic case 

For a completely disordered magnetic system, we have 

< >= (25) 
Therefore Xoi, Xn and Xoi vanish and we obtain for the energy 

E k /g = V(l ~ 27/0 - 2 - yr2(l- 7fc ) + [2-77(i + 2 7jfc )]2 (26) 

The dispersion relation Eq.(|26D is plotted in Fig.(fJ) for the physical value 
rj = 1/6. The band minimum is reached at points (±7r, ±7r). The center of 
the Brillonin zone T(0, 0) is a maximum. The bandwidth is given by 



W = E{0, 0) - E{tt, tt 
= 2g 



6+[l + ^P-l-^ 



(27) 
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for the physical value 77 = 1/6, we have W = 3.14(7. From Ref. we can get 
the relation between g and the t parameter of the t-J model, and we have 
t = (/(3 — 2rj)/A. We end then with a bandwith W ~ 4.711 We obtain a 
bandwidth of the order of t. We notice that even in this disordered case the 
bandwidth was reduced from 8t (for the t-J model with J = : the free 
motion) to 4.71i, which express the strong correlation in the system. We 
expect that in the ordered state the bandwidth will be reduced further. 

Now we want to examine the spin pattern around the hole. We started 
from a state without order and we want to calculate the local magnetization 
at site j, if the hole is at site i for the trial wave function ip^ with energy 
Ek- Algebric calculation leads to 

(^ k \riiS z ^) w k a k (r]e k - E k /g) 

T~ 7T\ r77T\ = —Wii COSKU — ?) 7i : . in (28) 

(^( fc ),^( fe )) 13 1 JJ e k a 2 k + 3[r)e k -E k /g] 2 1 ; 

where we have used Eqs.(P^D(P3|). Since the minimum of energy is at Q = 
(•7T, 7r), the hole will have the lowest energy Eq and then we get the following 
correlation function between the hole and surrounding spins 



For 1] = 1/6 we obtain 



E Q /g - 8 V 



6 + i[E Q /(7 - 877] 2 J 



Wij cos Q(i - j) 
(29) 



< riiS Z j >~ — - Wij cos Q(i — j) (30) 

Eq.(|5DD means that the hole polarizes the spins in its vicinity. The spin 
pattern around the hole is antiferromagnetic and the correlation decreases 
like when the distance between the hole and neighbor spin increases (see 
Fig.®). 

The trial state ip corresponds to a state with a total spin 1/2 and S^ ot = 
+ 1/2. The state corresponding to S% ot = —1/2 is just given by S^ ot ip. And 
the correlation function for a hole with spin down is related to the case of 
spin up by 



< ni S z 3 >^ = (S- ot 4>, ni S Z jS- tot ij) 

= - < n, t S z 3 > {T) (31) 
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This implies that a hole with spin down will polarize spins in its vicinity 
with opposite direction with respect to a hole with spin up. Naively we may 
say that when two holes with opposite spins are introduced in the system it 
will be preferable to have them close to each other in order to cancel their 
spin polarizations and obtain again a system without magnetic order which 
favours the motion of the holes. This scenario might be similar to the Cooper 
pairing. 



3.2 Locally ordered case 

Since experiments on copper oxides show that the spin susceptibility presents 
a sharp peak near the wave vector Q ||ll|| in the metallic phase, we shall take 
a lorentzian shape for the static susceptibility to describe the case of a locally 
ordered state. We write then 



Xq l + £ 2 [l + !(cos fe + cosg y )] 1 ' 

which may be correct only for short correlation lenght £ . C(£) is deter- 



mined by writing < SiSi >= 3/4 which gives 



c(0 = z 



iv? 1 



1 



g . + £ 2 [1 + |(cosg x + cosg„)]_ 
The inverse Fourier Transform of Eq. (132) gives 



iq(i-j) 



X\} ~ 3, 



1 v C(Qe 

N Y 1 + £ 2 [1 + \{cosq x + cosq y )] 

4 -^X Q (i-3)e- iQ{i - j) 



and Xq(^ — j) satisfies the equation 



(33) 



(34) 



J2 (~ A )u Xq(1 - j) + 7tXQ(i ~ j) = Si 
i s 



(35) 



with A the Laplacian on the square lattice. The asymptotic solution of this 
equation ( i.e. for | i — j \— > oo) is known and we have 
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where £l is given by 



(36) 



- = ch-\l + -) (37) 
For £ sufficiently large we have £l — £/v2- 

Now having the phenomenological form of the static susceptibility 
(Eq. (p2|)), we can compute Xoi, Xu an d Xoi for different correlation length. 
The table (1) contains some values of these quantities. We notice that by 
increasing £ we can have more and more ordered state. For example we al- 
most recover the values of Xoi, Xn an d Xo2 predicted by the linear spin wave 
approximation for £ = 5. 





Xol 


Xn 


Xo2 


1 


-0.105 


0.0278 


0.015 


2 


-0.199 


0.094 


0.058 


3 


-0.259 


0.150 


0.101 


4 


-0.299 


0.192 


0.138 


5 


-0.328 


0.223 


0.167 



Table 1 : Correlation functions Xoi, Xn an d Xo2 calculated 
from the phenomenological form of the static susceptibility 
Eq.(g|), for different £. 

Inserting the values of Xoi, Xn an d Xo2 in Eq.([H]) we get the dispersion 
plotted in Fig.(^|) for different values of £. We see from Fig.(||) that the 
spectrum changes drastically when the magnetic background changes from a 
disordered state to an ordered one. The minimum of Ek shifts from (±7r, ±7r) 
in the case of a disordered background (£ = 0) to (±7r/2, ±7r/2) for £ = 3. We 
see clearly that the bandwidth has been strongly reduced in comparison with 
the disordered case. It was reduced from W ~ 3. 14(7 ( f° r £ = 0) to W ~ 1.3g 
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( for £ = 3). In term of t we get a bandwidth W ~ 1.8t. But it is still of the 
order of t and not J. This is a discrepancy with previous calculations done in 
the framework of the t-J model for an AFM background |T], f|] which may 
be explained by the fact that no magnetization is present in our calculation 
since we considered our background as a singlet state. Other reason might be 
that with our variational approach we do not include the incoherent motion 
of the quasi-particle which is the consequence of retarded effects. 

4 Conclusion 

In this section we have computed the dispersion relation E{k) for a single 
hole in the framework of the spin-hole model using a variational method. 
The calculation was done for any rotational invariant magnetic background. 
And we have shown how the hole spectrum changes when the magnetic envi- 
ronment goes from a disordered background to an ordered one. Especially we 
have found that the minimum energy shifts from (±7r, ±7r) for the completely 
disordered case to (±7r/2, ±7r/2) in the case of ordered magnetic background. 
The bandwidth is strongly reduced by a factor 1/3, when the magnetic back- 
ground changes from a disordered to an ordered state. An other interesting 
result that we have found is the spin polarization around the hole. Namely 
when a hole moves in a completely disordered backgound it will polarize spins 
in its vicinity. This polarization is decreasing with the distance as a power 
law and has an AFM ordering. We expect that two holes with opposite spins 
will tend to pair in order to remove their respective spin polarizations. This 
question deserves further studies and it is now under investigation. 
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Figure 1: Dispersion relation for a single hole in disor- 
dered magnetic background( with M = (n, ir) , Y = (0, 0) 
and X = (it, 0) ). is in unit of g. 



Figure 2: Polarization of spins surrounding the hole. 
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Figure 3: Dispersion relation for a single hole in differ- 
ent magnetic backgrounds : from a completey disordered 
state — 0) to an ordered state £ = 5. Ek is in unit of 
g ( with M = (tt, tt) , T = (0, 0) and X = (tt, 0) ). 
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